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Abstract: We compute the short distance coefficients for the NRQCD factorization for- 
mulas of the meson-to-vacuum matrix elements for the axial vector and vector parts of the 
charged weak current in the 5-wave spin-singlet and -triplet be mesons, respectively. The 
computation is carried out to order including relativistic corrections of all orders in q^", 
where q is the relative momentum of the b and c in the meson rest frame. The relativistic 
corrections at order Og are new. The results reveal that the relativistic corrections to the 
leptonic decay rate of the Be meson at order or less converge rapidly, which shows a 
strong contrast to the uncomfortably large corrections of order a^lgl*^. The short distance 
coefficients listed in this paper can be employed to compute the resummation of relativistic 
corrections to the phenomenological measurables that involve and B* production and 
decay. 
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1. Introduction 

Among various quarkonia including charmonia and bottomonia, the bound state of a be 
pair is a distinct heavy quarkonium system composed of two different heavy quark flavors. 
As the spin-singlet 5- wave bound state of the be pair, the Be meson was discovered by the 
CDF Collaboration at the Fermilab Tevatron through the decay mode Be — s- J / ip + + i^i 
im and the 5- wave spin-triplet state B* has not been observed, yet^. The be meson offers a 
unique laboratory for the nonrelativistic quantum chromodynamics (NRQCD) factorization 
framework Q because the typical heavy quark (Q) velocity vq in the be meson lies between 

^B^ is the bound state of a be pair and is that of cb. Throughout this paper, Be denotes B^ . 
However, our analysis can be equally applied to the charge-conjugate state B~ . 



- 1 - 



that (w^ ~ 0.3) in the charmonium system and that (u^ ~ 0.1) in the bottomonium system. 
In addition, unhke the spin-triplet 5- wave QQ mesons, J/V' and T that can decay into 
lepton pairs through the electromagnetic current, the decay of the Be meson can proceed 
only through the charged weak current and various dynamics play their roles in the decay 
of the he meson |^ § ||, 0, §]• As of now, except for the mass, niB^ = 6.277 ± 0.006 
GeV, and the life time, tb^ = (0.45 it 0.04) x 10~^^ s, for the spin-singlet S'-wave state 
Be, little is known about various properties of the be bound states experimentally Q. In 
near future, one may probe many unknown properties of the bound states in detail as the 
CERN Large Hadron Collider accumulates orders of magnitude larger number of events 
than the currently available data for the be bound states. Therefore, it is desirable and 
necessary to achieve better accuracies in theoretical predictions. 



Earlier theoretical studies on the be bound states cover the spectroscopy [11, 12, 13 



14, 15 1 and the production mechanism at colliders ||T^, 17, |l^, In order to achieve 
better accuracies in the predictions for the 6c-meson production and decay, it is necessary 
to know accurate values for the decay constants. The NRQCD factorization formula is 
useful in making a systematic series expansion of the decay constant for the be meson 
system in powers of vq. The NRQCD factorization theorems have been proved for the 
electromagnetic and light hadronic decays of heavy quarkonia [|| and for a few exclusive 
production processes of heavy quarkonia |2^, 23|. The NRQCD factorization formulas 
for the hadronic part of the charged weak current that involves the decays of Be and B* 
are similar to those for the electromagnetic decays of QQ mesons. Braaten and Fleming 
computed the one-loop QCD corrections to the Be decay constant in the static limit vq = 
and the relativistic corrections of relative order a^Vg [24| by computing the short distance 
coefficients of the NRQCD factorization formula for the decay constant. Based on the same 
strategy, Hwang and Kim calculated the B* counterparts [25|. Two-loop QCD correction 
to the short distance coefficient for axial vector current involving Be decay constant was 
calculated by Onishchenko and Veretin which shows uncomfortably large correction 
like those for the leptonic decays of the spin-triplet S'-wave quarkonia J/ip and T |2^, 
In fact, according to the velocity-scaling rules of NRQCD [Q], the corrections of relative 
orders Vq and asVg should be equally important as that of relative order a'^Vq. In addition, 
the large separation between nii, and nic, where mq is the mass of the heavy quark Q = 
b or c, in the be system gives rise to factors of log{nib/mc) in the one-loop corrections 
1 24, 25, 26 1, which may potentially deteriorate the convergence of the power expansion in 
Vq. Therefore, it is worthwhile to check if such large corrections indeed arise when one 
includes the relativistic effects. 

In this paper, we compute the relativistic corrections to the NRQCD factorization 
formulas of the meson-to-vacuum matrix elements for the axial vector and vector parts of 
the charged weak current in the S'-wave spin-singlet and -triplet be mesons, respectively. 
The calculation is carried out at first order of the strong coupling including relativistic 
corrections to all orders in vq. The short distance coefficients for the NRQCD factorization 
formula are usually obtained after subtracting the infrared (IR)-sensitive contributions 
of the NRQCD correction from the full QCD corrections by perturbative matching. It 
requires laborious bookkeeping of the Feynman rules in the NRQCD perturbation theory 
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which grows tremendously extensive as the order in Vq" increases p9| , |30| ] . Instead, we use 
a new method introduced recently in to compute the relativistic corrections to the 
short distance coefficients at order covering all orders in v'^. The method integrates 
out the temporal component of the loop momentum by contour integration in the one-loop 
corrections to the full QCD amplitude and, then, expands the integrands in powers of 
the external momenta divided by mg and the spatial components of the loop momentum 
divided by mq. In comparison with, so called, the method of region given in p^ , this 
new method is particularly useful in computing relativistic corrections of higher orders in 
vq and even makes it possible to find the closed form of the expression that includes the 
relativistic corrections resummed to all orders in fg" at one loop. 

This paper is organized as follows. In section |2|, we discuss the perturbative match- 
ing of NRQCD onto QCD at one loop. Kinematics of the problem and the definitions 
of the variables that are useful in computing the short distance coefficients are given in 
section ^. Section Q contains the strategy and detailed formulas to compute the short dis- 
tance coefficients. We compute the QCD one-loop corrections in section |5| followed by the 
NRQCD corrections in section |^. Our final results for the short distance coefficients are 
listed in section ^ and we summarize in section ^. In appendices, we provide the formulas 
for the tensor-integral reduction and list the values for the loop integrals that appear in 
the one-loop correction to the QCD and NRQCD amplitudes. 



2. Perturbative matching to all orders in v 

We define the hadronic parts of the weak decay amplitudes iA'i^ and iA'^* for Be and B* 
as the meson-to- vacuum matrix elements of the axial vector current b^^'j^c and the vector 
current 67^0, respectively as^: 

i^5,B. = (0| h^l^clB,) = ifB^K'^, (2.1a) 
= (0| h^clB*,) = ifB-imB-iS^, (2.1b) 

where h and c are the Dirac field operators for the bottom quark and charm quark, respec- 
tively. The amplitudes (^]^) are scaled by the leptonic decay constants fu for H = Be 
and B*. Here, K is the meson momentum, e is the polarization vector of the B* meson, 
ms, = 6.277 (6) GeV @ and uib* = 6.330 (7)(2)(6) GeV [|| are the masses of the B^ 
and B* mesons. The quarkonium state \H) in ( |2.1| ) for H = B^ and B* is normalized 
relativistically: {H{K')\H{K)) = 2 K° {27r)^ 6^^'> {K' - K). In the meson rest frame, the 
matrix elements in ( p.ip become simple: Because K = {mH,0) in this frame, only the 
0-th component survives in the matrix element ( |2.1a| ) for the B^.- Due to the transverse 
condition e ■ K = for the B*, e = (0,e) in this frame and, therefore, only the spatial 
components are nonvanishing in the matrix element ( ^.Ib ) for the B*. 



^See, for example, 
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According to the NRQCD factorization Q, we can write the nonvanishing components 
of and iA^* in the rest frame of the be bound states as 

zAIb^ = y^2^^Y.Pn{0\Onm, (2.2a) 

n 

iA%, = v/2^j;K(0|O;|i?:), (2.2b) 

n 

where P„ and Vn are the short distance coefficients and On and are the NRQCD 
operators. The operator matrix elements in ( |2.2| ) are regularized dimensionally in d = 4— 2e 
space-time dimensions. The overall factor \J2mi{ for H = Be and B* in ( |2.2D has been 
taken out because the state \H) in ( p.2p is normalized nonrelativistically: {H{K')\H{K)) = 
{2t^)^ 5^'^\K' — K), while the amplitudes on the left sides of (2.2) have the relativistic 
normalization for the quarkonium H like those in ( |2.1| ). 



The main purpose of this paper is to compute the short distance coefficients Pn and 



Vn in (2^) that correspond to be color-singlet operators at order as- These coefficients can 



be determined by the matching equations 

iKbc^=Y^Pn{0\On\bei), (2.3a) 

n 

i4c,=Y.^n{0mbei), (2.3b) 

n 

which is the statement of NRQCD factorization for the perturbative color-singlet be state. 
Here, 6ci denotes the color-singlet be pair whose invariant mass is the same as the meson 
mass. Throughout this paper, we suppress the factor ^/T^ that comes from the implicit 
color trace in iA^ t and iAl , where A'^f = 3 is the number of colors. Note that the 



5,6ci be 



coefficients P„ and Vn in (O) are identical to those in ( p. 21 ) because the short distance 
coefficients must not depend on the long distance nature of the heavy quarkonium state. 
While the amplitudes (|2.2| ) contain nonperturbative quantities, the amplitudes (^) are 
calculable perturbatively. However, it is possible that the amplitudes (|2.3|) acquire singu- 
larities in the IR or ultraviolet (UV) regions at order Og or higher. These divergences are 
to be regularized dimensionally. The matching equations that contain the terms upto order 
as are 



n n 



where the superscripts (0) and (1) indicate the order in Ug. In the first sum of each line in 
( |2.4| ), only color-singlet be operators contribute, while in the next sum additional operators 
may enter once they mix with color-singlet be operators under one-loop QCD corrections. 
Through order a<j, the NRQCD amplitudes are defined by 

MsgjNRQCD = Y.Pi'H0\On\be,)^'\ (2.5a) 

n 

M£Vrqcd = E K^°^0|O;|6ci)(^-), (2.5b) 
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where j = or 1. At order a°, the NRQCD matrix elements (0|O„|6ci)(°) and (0|C'J,|6ci)(°) 
are finite and the short distance coefficients -P^^^ and Vn^^ can be determined from the 
identities 

M°;£JnRQCD = ^Af^^^ = Y: Pi'H0\O^.\-bc,f\ (2.6a) 

n 

M-,i?]NRQCD = = E Vi'^Homlcfl (2.6b) 



Each expansion of (2^) is a power series in q/mq, where q is half the relative three- 



momentum of the b and c in the center-of-momentum (CM) frame of the be pair. However, 
at order Og, the short distance coefficients P^^^ and Vn^^ must be determined after sub- 
tracting the long distance contributions that are contained in the order-a^ matrix elements 
(0|O„|6ci)^"^^ and (010^1601)^"^^ which include the contributions from the potential, soft, 
and ultrasoft regions: 

^A'J^yl^ - M°£]nrqcd = E Pi'H0\On\bc,)^'\ (2.7a) 

n 

- Mg^RQCD = E Vji'Hombc^f\ (2.7b) 

n 

In this way, one can determine the order-a^ short distance coefficients pji'^ and Vn^\ which 
are free of IR sensitivity. 

The computation of [iA^J"^^^ ]nrqcd and [i^g^^^'']NRQCD is very complicated because 
it involves operators and interactions of all orders in vq. Fortunately, the authors of 
1 31] recently introduced a way to compute the NRQCD amplitudes [^^^''^^^Jmrqcd and 

[i^^^^^^'' Jnrqcd directly from the full QCD counterparts iA^^^^^^ and iA^^^ by expanding the 
integrands in powers of the momentum divided by mq. Before the expansion, temporal 
component of the loop momentum has been integrated out, using contour integration, in 
order to avoid the generation of ill defined pinch singularities that may develop if one ex- 
pands the heavy quark propagators too early. The integrands for the integration over the 
remaining spatial components of the loop momentum are then expanded in powers of the 
external momenta divided by mq and the spatial components of the loop momenta divided 
by mq. Divergent integrals over the spatial components of the loop momenta are regu- 
larized dimensionally while scaleless power-divergent integrals are dropped in accordance 
with the dimensional regularization scheme. As the last step, remaining UV divergences are 
renormalized according to the modified minimal subtraction (MS) scheme. In this work, we 
apply this method to compute the NRQCD amplitudes [i^^^j^^ Jnrqcd and [i»4^|,^^^]NRQCD 

and determine the short distance coefficients Pn^^ and V^^^ from ^ 
3. Kinematics and notations 

In this section, we define notations for the kinematics of the problem. We take pc and pi 
to be the momenta of the incoming heavy quark c and heavy antiquark 6, respectively. 
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which are on their mass shells: = and 



mt 



They are expressed as linear 



combinations of half the total momentum p = \{pc + h) ~ \^ ^'^'^ hsli their relative 
momentum q = \{pc — Pi)'- 



Pc = p + q, 
Pb = P-(l- 

In the CM frame of the be pair, the momenta are given by 

Pc = {Ec,q), 

Pb = iEb,-q), 
p = [^{E, + Eh),0], 
q = [\{E,-Eb),ql 

where Ec = {ml + q2)i/2 ^nd Eb = {ml + 9^)1/2^ Note that 

Ap-q={pc+ Pb) ■ {Pc - Pb) = ml - ml 0, 



(3.1a) 
(3.1b) 



(3.2a) 
(3.2b) 
(3.2c) 
(3.2d) 



(3.3) 



unlike the case of the QQ pair considered in |31|. 

For later use, it is convenient to define parameters 5 which is the magnitude of the 
three-momentum of 6 or c and Cc which is the energy of the charm quark scaled by the 
invariant mass y^4p2 = + Eb of the be pair in the CM frame: 



Er 



{mb + mc 
,2 ^2 



{mb - mcY^ 
4p2 



1 



/4p2 2 V 4p2 

Note that 5/ec = \q\/Ec. The following relations are also useful: 

Pl-Pc + m^mc = {EbEc + m^mc) + q^ = 2mbmc + 0(q^ 

{mb + mc)^ 



Pb-Pc 



mbmc = {EbEc - mbmc) + q 



2mbmc 



q' + 0(9^ 



(3.4a) 
(3.4b) 

(3.5a) 
(3.5b) 



In the derivation of the full QCD amplitudes, we use the Dirac spinors for the c and b 
with the nonrelativistic normalization. In the CM frame of the be pair, they are 



Uc{Pc) 
Vc{Pb) 



Mb 



{Ec + mc)ic 
q-cr(.c J 

-q ■ o-'nb 

{Eb + mb)'nb 



(3.6a) 
(3.6b) 



where Mq = [2EQ{EQ + mQ)]~^^'^ for Q = c and b, and r]b are Pauli spinors for the c and 
b, respectively. The spinors in ( |3.6D are convenient in making nonrelativistic expansions. 
The threshold expansion method in 1 34 ] and its dimensionally regularized version in |^5| also 
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use this form except that the relativistic normahzation is used. To extract the spin-singlet 
and -triplet states from a full QCD amplitude for the be pair, one can also make use of 
the spin-projection operators for those states. In []36| , |23[| , for example, one can find the 
spin-projection operators for the spin-singlet and -triplet states of a heavy quark-antiquark 
pair with different flavors. 

We use the representation for the Dirac matrices introduced in [51]: 



where 1 is the identity matrix. In ( |3.7D 7' and the Pauli matrix are defined for z = 1, 2, 
• • • , d — 1. The requirement of the Clifford algebra for the Dirac matrices in d space-time 
dimensions 

{7^7"} = 25^^1, (3.8) 
for n, 1/ = 0, 1, 2, ■ ■ ■ , d — 1 forces the anticommutation relations for the Pauli matrices, 

{a\a^} = 25'n, (3.9) 

for i, j = 1, 2, • • • , d — 1. In our computation of the spin-singlet case, we encounter the 
loop correction to the axial vector current. We carry out the Dirac algebra by making use 
of the naive dimensional regularization, in which 75 anticommutes with for any indices 
/X in d dimensions. This prescription is self consistent for the case considered in this paper 



1 37]. As commented in section 3 of [38|, the matrix representation of the Dirac 75 that is 



consistent with the choice (3^) is then 

75 = C : 1 • (3.10) 




which guarantees {7'^, 75} = 0. With the matrix representations for the spinors in (|3.6| ) and 
with the set of Dirac matrices (|3.7| ) and ( |3.10| ) , we can carry out the calculation consistent 
with naive dimensional regularization. 

For d — 1 = 3, the Pauli matrices satisfy the commutation relations 

[a\a^ = 2ie'^''a''. (3.11) 

However, for the spatial dimensions greater than 3 the totally antisymmetric combina- 
tion of three Pauli matrices {[a^,(T^],a^}, which may arise in the threshold expansion of 
the products of three or more Dirac matrices of different spatial indices [^], is linearly 
independent of both 1 and o"^ for £ = 1, 2, • • • , d — 1. And the reduction 

{[a\a^],a''} = Aie'^''l (3.12) 



is allowed only at 3 spatial dimensions ^5[. Therefore, unless divergent contributions 
disappear, we do not use the reduction ( |3.12 ). 
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4. Formulas for short distance coefficients 



Let us first classify the operators tliat appear in the matching conditions ( |2.6[ ) and (p.Tj). 
For the spin-singlet 5-wave case, only a single type of operators On contributes and for the 
spin-triplet S'-wave case, there are two kinds of operators and 



An 



Bn 



xl{- 



\2n „i 



(4.1a) 
(4.1b) 
(4.1c) 



where is the Pauli spinor field that annihilates the charm quark and is the Pauli 



spinor field that annihilates the antibottom quark. The operators in (4J.) contain ordinary 
derivatives, rather than covariant derivatives so that they are not gauge invariant. We 
evaluate their matrix elements in the Coulomb gauge, in which inclusion of the he operators 
involving the gauge fields brings in corrections of relative order Vq |23|. While the operators 
On and 0\^ have only the 5-wave contributions, the operator also contains the D- 
wave contribution as well as the 5-wave one. The operator can be decomposed into 
a linear combination of 0\^ and the D-wave operator 0\)^: 



O 



Bn 



d-l 



0\n + 



Dni 



(4.2) 



where C})„ is defined by 



^hn 



xli-h V) 



2n-2 



v)(- 



iv) 



d-l 



(-1 v)V^ 



(4.3) 



In the basis of operators 0\^ and for the spin-triplet case, the matching conditions 
dU) and (|]7|) become 



4? = E«^^^^<0|^^jSci)W + j:6(°)(0|O^„r6ci 



,(0) 



(4.4a) 



bci 



iX 



bci 



NRQCD 



j;aW(0|O^J6ci)W + J^6«(0|O),J6ci)(°), (4.4b) 



where and 6„ are the short distance coefficients corresponding to the operators 0\^ 
and O^^, respectively. A similar equation holds in the basis and O^^^, where the 
associated short distance coefficients are 



a-n + 



d-l 



bn, 



n — bn. 



(4.5a) 
(4.5b) 



where Sn and Dn are the 5-wave and D-wave components of the short distance coefficient 
Vn, respectively. For more details, see [31|. 
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The be matrix elements for the spin-singlet case in (|2.6a ) and ( p. 7a ) and those for the 
spin-triplet case in ( p.6b ) and (|2.7b| ) are calculable perturbatively as 



(0|On|6ci)W 
(0|Ok|feci)(°) 
(0|Ol,J6ci)(0) 



(4.6a) 
(4.6b) 
(4.6c) 



where and ijb are two-component spinors for the charm quark and the antibottom quark, 
respectively. In order to maintain consistency with our calculations in full QCD, we have 
taken the be states to have the nonrelativistic normalization and we have suppressed the 
factor x/Nr that comes from the color trace. 



The most general Lorentz covariant forms of iA'^r and iAr are 



5,bci bci 

for the axial vector and vector currents, respectively. Here, 



5, feci 

iAf^ 



G = yz^(l + A), 



(4.7a) 
(4.7b) 

(4.8a) 
(4.8b) 



where Zq for Q = h ox e\s the wavefunction renormalization factor of the heavy quark Q 
at order |24]: 



1 



27r 



4euv 



1 3 , 47r/x2e-TE 

+ + - fog + 1 

2eiR 4 mjj 



(4.9) 



where Cp = {N^ — l)/{2Nc) = 4/3, ^ is the renormalization constant, 7e is the Euler- 
Mascheroni constant, and the subscripts on 1/e indicate the origins of the singularities. A5 
and A in ( |4.8| ) are the multiplicative corrections to the axial vector and vector vertices, 
respectively. Note that the terms proportional to p'^ survive in ( |4.7D because the weak 
currents are not conserved while those terms vanish in the electromagnetic current which 



is conserved. At order a^, the only nonvanishing contributions in (|4.7| ) are G 



(0) 



and all the other contributions are absent: H, 



(0) 



i?(0) = g{0) = 0. The leading 



nonvanishing contributions to H^, Q^, H, and Q appear from order ag- 

Similarly, the nonvanishing components of the NRQCD counterparts to iA'^j^^^ and 
iA^^ in the CM frame of the be pair are 



A"- 



NRQCD 



Vb{p-b){G 5, NRQCDl° + -H's.NRQCDg" + Q5,NRQCDP°)75 Itc(Pc) , (4.10a) 



A 



feci 



^fe(P5)(G'NRQCD7' + Htimci)q')uciPc 



NRQCD 

where we have used = in the CM frame of the be pair and 

G5,NRQCD 

Gnrqcd 



Zb]^RQCDW ■^c]nRQCd(1 + As^nRQCd), 
■^fe]NRQCD[\/^]NRQCD(l + Anrqcd)- 



(4.10b) 



(4.11a) 
(4.11b) 
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The heavy quark wavefunction renormaUzation m NRQCD, [-^ Zq ]nrqcd for Q = b or c 
at order a^, is given by pl| 



I NRQCD 



1 + 



UsCf 



1 



2vr V euv eiR 



(4.12) 



By making use of the Dirac spinors in (|3.6D and the Dirac matrices in (3/7) and ( 3.10 ), we 
find that the expressions in (13) and ( p^ are reduced into 



Vb{p-bh^l5Uc{Pc) = MbMc[{E^ + mb){Ec + rUc) - q^] rjUc, (4.13a) 

VbiPbh^UciPc) = -MbMc[{Eb + mb){Ec + m^) + q^] r^ic, (4.13b) 

Vb{Pb)l'uc{pc) = MbMc{[{Eb + mb){E, + me) + q^] vI^t'Cc - 2q' vU ■ <t^, (4.13c) 

Vb{pi)uc{pc) = -MbMc{Ei, + mb + Ec + nic) rjlq ■ a^c, (4.13d) 



where the expressions in ( 4.13| ) are valid to all orders in q and we have used the identities 



q', 



{q . af 

(q ■ CT)a\q ■ a) = 2q\q ■ a) - a'q^ 



(4.14a) 
(4.14b) 



that derive from ( |3.9| ) in d — 1 spatial dimensions. Note that we do not encounter the 
products of Dirac matrices which involve products of three Pauli matrices of different 
indices that bring in the contribution in ( 3.12| ). Therefore, the threshold expansion ( 4.13| ) 



is free of ambiguities in d — 1 spatial dimensions. Substituting ( 4.13 ) into (|4.7| ), we find 
that 

'^Im = ^b-^cjcs [{Eb + mb){Ec + rric) - q^] 

-{H^q"^ + Q^p^) [{Eb + mb){E, + m,) + q^] jr/^^^c (4.15a) 
= NbMc\G[{Eb + mb){E, + me) + q^]r]yic 

- [2G + {Eb + mb + E, + m,)H] q'rjl{o- ■ q^^ (4.15b) 

Similarly, the NRQCD amplitudes [^.4^ ^^Jnrqcd and [«.4^^Jnrqcd in (jl.lOl ) are the same 



as iA^T. and iAl in (4.15) except that the coefficients F and F^ are replaced with 



5,f)ci feci 

i^NRQCD and Fs^NRQCD, respectively, for F = G, H, and Q. 

According to the matching conditions in ( |2.6D , [^.A^^^^^Jnrqcd and [^-4.^111^ Jnrqcd are 
identical to the order-a^ full QCD counterparts in ( 4.15 ). In order to obtain the short dis- 



tance coefficients at order a^, we need to expand [i^g^^]NRQCD 



^-^5^1 ^"^d M^^ JNRQCD 

.^^^ as linear combinations of the perturbative NRQCD matrix elements in ( ^4.61) . By 
making use of these order-a^ perturbative matrix elements and the matching conditions in 



li{0)i 



iX 



m 
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^^), ( [4.4a ), and (|4.5| ), we can obtain the short distance coefficients at order a^: 



1 / f) \ ^ 



rrir 



q{0) 



<j2=0 



(4.16a) 

(4.16b) 
(4.16c) 



where we have used the fact that = G^") = 1 and = = = Q(o) = and 
the short distance coefficients Un^ and bn^ are given by 



n! 



5 \ " 



(n — 1)! \9q 



9 



ra-l 



<j2=0 



£?2=0 



(4.17a) 
(4.17b) 



At order al, the quantities A5, A, and Zq in ( [4.8[ ) contain divergences. The multi- 
plicative vertex correction factors A5 and A have logarithmic divergences in the UV and IR 
regions. They also contain Coulomb divergence, which is not analytic in the limit |q| — )• 0. 
The wavefunction renormalization constant Zq has logarithmic divergences in the UV and 
IR regions. However, because of the usual cancellation between the vertex and fermion 
wavefunction renormalizations, G^^ and G^^^ are free of UV divergences and contain only 
IR divergences. The quantities H^^\ H^^\ and Q^^\ which contribute only from this 
order, may have only Coulomb divergences in the limit \q\ — )• 0. Therefore, in the full 



QCD amplitudes iA^J'j^^^^ and , UV divergences cancel and the amplitudes may have 

singularities only in the IR region, which are either Coulomb or logarithmic divergences. 
Because NRQCD reproduces full QCD in the IR region, the IR divergences in Gj^^j^qqj-,, 



(1) 



5,NRQCD' V5,NRQCD' '-'NRQCD' ^^^"-^ -'-'NRQCD 

H^^^ , respectively, to make the following quantities IR finite: 



G 



(1) 



and H, 



(1) 



cancel those in G^, m^' , , G^^\ and 



AG 



(1) 
5 

(1) 



G 



ag(i) 



(1) 



G 



(1) 

5,NRQCD' 



Q 



5,NRQCD' 
(1) 

5,NRQCD' 



G(^) - G 



(1) 

NRQCD' 



(1) 

NRQCD- 



(4.18a) 
(4.18b) 
(4.18c) 
(4.18d) 
(4.18e) 



Therefore, the right sides of ( ^.7] ) which are determined by the expressions in ( 4.18D are 



free of IR sensitivities. Now we can obtain the short distance coefficients at order a], by 
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making use of the matching conditions in (^.T]), ( [4.4b ), and (4^) as 
Pi'^ = ^ AA,AAe| AG^') [{Eh + m6)(i5;e + me) - g^] 



q(i) 



5(1) 



where the short distance coefficients ^ and bn^ are given by 



,(1) 



5(1) 



1 



5 



1 



d 



(4.19a) 

(4.19b) 
(4.19c) 

(4.20a) 



92=0 



(n - 1)! \dq^ 



n-l 



(4.20b) 



Although we have completely removed the IR singularities in the short distance coefficients 
( |4.19 ) and ( 4.20| ), they have logarithmic UV divergences that are originated from the one- 



loop NRQCD matrix elements in Gg^^j^Q^^j-, and G[jp^Q(-;j-, . As stated before, the quantities 
AhP'\ Aq1^\ and AH^^^ are free of UV divergences as well as IR divergences. We 



renormalize AG^^^ AG^^^, and the short distance coefficients in ( [4.19 ) and ( [4.20| ) according 
to the MS scheme to find that 



[f (1)1—= f (1)1 

IJn J MS 1aF{i 



(1)1 



(1)5 



(4.21) 



where F = G5 or G and f = P, a, b, S, and D. In deriving the expression for [/i^^] in 



( 4.21 ), we have used the fact that, in minimal subtraction, one removes the 1/e pole times 



the order-Og d-dimensional matrix element. Hence, a term proportional to {d — 1) 



-1^-1 ; 



IS 



subtracted in (4.19b) in carrying out the renormalization. 



5. QCD corrections 

In this section, we compute the one-loop QCD corrections to the axial vector and vector 



parts of the charged weak current iA^^^^^^ and iA^'^J^' , respectively. The order-a^ QCD 
corrections are composed of the vertex corrections and the wavefunction renormalization 
contributions. In the Feynman gauge, the vertex correction contributions to iA^^^^^^ and 

lA^^ are given by 



0(1) 



AM 



-iglCp 



k D0D1D2 
DoDiD2 



(5.1a) 
(5.1b) 
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where g1 = AttUs is the strong couphng and the symbol and the denominator factors 
DiS are defined by 



Ik 



Do = k'^ + ie, (5.2b) 
Di = -2k-pi + i£, (5.2c) 
D2 = + 2k-pc + ie. (5.2d) 

Here, /x is the renormahzation scale and we have taken the gluon momentum as the loop 
momentum k. By making use of the anticommutation relation of 75 and 7^^ we rearrange 
each term of the numerator in ( ^.l] ) as a linear combination of terms '^a^n"l5 or 7ari7", 
where Fj is a product of three or less Dirac matrices. Applying the anticommutation 
relation (|3.8| ) summed over the d-dimensional index a and making use of the on-shell 
conditions 

i)cu{pc) = mcu{pc), (5.3a) 
k^iPb) = -^bvipi), (5.3b) 

we can reduce the expressions in (0) as 

AM ■ 2n f v{p-bW^ik,P,q,mh,-mc)'y5u{pc) , . 

^5 = -^SsCf , (5.4a) 

= -^glC, ^ ^(H)rnfc^P,^,rn,,m.)^(p.)^ ^^^^^^ 

where T>^{k,p,q,mh,mc) is defined by 

r>'{k,p,q,mb,mc) = [{d - 2)k'^ - 2{4p'^ - ml - ml) + 8k ■ q]-/'' 

+2{mbi^^ + mc^-i^') + 2(2 - d)k^'^ - Sgf"^. (5.5) 

By making use of the standard reduction methods for the tensor loop integrals, we can 
express all of the loop-momentum dependence in terms of pi and pc, which are linear 
combinations of p and q. Because p ■ q ^ 0, the reduction formulas for the tensor integrals 
of the S'-wave bci decay are slightly more complicated than those for the spin-triplet S- 
wave QQi decay in [^] where p ■ q = 0. In appendix ^ we list the formulas for the tensor 
reduction. Once we apply the equations of motion in (|5.3|), F'^'s in ( |5.4D are reduced into 
a linear combination of 7^, p^l, and q^l as 

= -iglCFvip-b)[M-mch'' + A2i-mc)p'' + A3{-mc)qf']j5u{pc), (5.6a) 
A^ = -iglCpviPi) [Ai{m,)j'' + A2{m,)p^ + ^3(^^)9^] n(pe). (5.6b) 
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As mentioned in [p5| ], the Lorentz scalar coefficients Aj's in ( ^.6a ) and ( ^.6b ) are the same 
except for the replacement Ai{—mc) -H- Ai{mc) and Ai^mcYs are defined by 



Ai{mc) 
Mime) 



{d - 2)(Ji - 2J4) - % • pcJ2 + 4J3 + 
2 



(mfe - rricYJz + (mf - m^) J5 



P5 • Pc + Jn-feWic 



(mfe - mc) J3 + (mfc + m-c) J5 _^ (2 - d)JT{mc) 



Pb-Pc + mbUic 



2p^q 



.2r,2 



{nib + "T-c)'/3 + ("T-b - "T-c) J5 _^ (2 ~ d)J6{mc) 



Pb-Pc 



mbruc 



2q^ 



(5.7a) 
(5.7b) 
(5.7c) 



Here, the Jj's are scalar integrals which are defined and evaluated in appendix ^ Except 
for the two scalar integrals Jg and J7, the other scalar integrals Jj for i < 5 are even 
functions of rric and, therefore, we suppress the arguments of Jj for i < 5. The expressions 
in (|5.6|) are new, in which the relativistic corrections to all orders in are included. 



According to ( 4.7a| ) and (|4.8| ), the vertex corrections can be parametrized as A^ = 
v{pi){K^'~^'^ + H^p>^ + Q^q^^)^^u{pc) and A'^ = v{pi){K'~^>^ + Hp^ + Qq^)u{pc). Therefore, the 
multiplicative factors A5 and A for the vertex corrections can be determined by ^i(=f?n,c), 
which are the coefficients of in ( [5.6[ ). In similar ways, and H are determined by 
A2{^mc) and Qs and Q are determined by ^3(=Fmc), respectively. Substituting the values 
for the scalar integrals Jj's evaluated in appendix ^ into Ai{^mc) in ( |5.7D and then substi- 
tuting Aii^rric) into (^.61), we obtain the multiplicative factors for the vertex corrections 
as A5 = —ig1CFAi{—mc) and A = —iglC f Ai{mc) ■ The results are 



A, 



asCp] 1 1, (4Vg-7E)2 

— ^ < + o log o-^^ + 

4vr euv 2 



Pb-Pc 



1 1 , (47r/x2e^TE)2 
— + - log ^ 



2p2 VeiR 2 



+ <5M 6 + 



2(mb + mc 



rribmc 



+2pi-pcL2{6, Cc)] log— I + 



mt 



2p^ 



A = A 



5 Irrtc 



p2 



9 ■ 4 2 2 

^ 25 °^ [(p5 • pcY - m^m^i 



(5.8a) 
(5.8b) 



where A^lm^^-mc represents the expression that can be obtained from A5 after replacing 
rric with —rric and the variables 6 and Cc are defined in (^^). The functions Li((5, Cc), 
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L2{S,ec), and K{5,ec) are defined by 



Li{5, Be) 



12(5, Be) 



Ki6,ec) 



J_j (ec + S){l-ec + 5) 
26 °^ (ec -6){l-ec- 6) 

{rrib + mc)^ (m;, + rricf {ml - AmcUib + ml) 2 , 4^ 
mbmc <^-^^.i-^^.-. 

Ij (ee + (^) (1 - e, - (5) 
26 °^ (ee - <5) (1 - e, + 5) 



{mb + mc){mb — mcY 



mbmc 



6mfm^ 



1 

45 



Sp 



25 



Sp 



-25 



+ Sp 



25 



(5.9a) 



(5.9b) 



ec + 6 J \ec — 6 J \1 — Cc + 6 
{mb + mc)^ (?Tifc + rric)^ (m^ - l^m^mb + m^) 



Sp 



-25 



q' + 0{q-). 



(5.9c) 



We have listed the first two leading terms in the nonrelativistic expansions of these func- 
tions. According to (|5.9|) , the functions Li(5, Cc), ^2(5, ec), and K{6,ec) are of order 1 as 
\q\ — )• 0. The function Sp(2;) in ( ^.9c ) is the Spence function, which is defined by 



Sp(x) 



log(l - t) 



dt. 



(5.10) 



Now we evaluate G5 and G. Substituting A5 and A in (|5.8| ) and the heavy quark 



wavefunction renormalization constant Zq for Q = b and c in (4^) into (|4.8| ), we obtain 



1 + 



47r 



^(^Li(5,eJ--)-2 



1 1 , (47r//2e-TE)2 
eiR 2 



miml 



+ 



6 + 



2(mf, + mc) 



2 1 



Pb-Pc- mbmc 
+2pi-pcL2{6,ec)] log^ + 



-^l('5, Cc) - y 



52 

m\ 



^K{6,ec 



1 



2p^ 



vr 



log 



2 2 



4 2 2 

p mj^m^ 



5 25 [(Pb-Pc)^ 



(5.11a) 
(5.11b) 



As we have expected, G5 and G are free of UV divergences while they have Coulomb and 
logarithmic divergences in the IR region. 

Next we substitute the values for the scalar integrals Jj's evaluated in appendix ^ 
into A2{^mc) and A3(=Fmc) in (^). Then we can determine = —iglGF-A2{—mc)j 
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-ig1CFA?j{—mc), H = —ig'^CFA2{mc), and Q = —ig1CFA2,{mc). The results are 



Q5 



H 

Q 



rrir 



4tt Pi-Pc + marine 

UsCf 



Li{5, Cc) - y 



+ r_ log — 
ml 



-{mb + mc 



Pb-Pc- rabmc 



Li{6, Cc 



ITT 

T 



+ r+ log ■ 



mn 



Li{6, ec 



ITT 

T 



+ r_ log 



p [ 



1 + r+ log — I 



-rUc ) 



(5.12a) 



(5.12b) 

(5.12c) 
(5.12d) 



where r± is defined by 



r± 



■mb =F iric 
4:{mb ± mc) 



(mb ± mc 
4p2 



(5.13) 



The variable r_ is of order 1 and the variable r+ is of order q^, which are both finite 
in the limit |g| — t- 0. The real parts of the quantities in ( ^.12 ) are finite. According to 
( |3.5| ) and (|3.4[ ), luiH^ is of order |g| and, therefore, finite and the leading contribution to 
ImQs^NRQCD is of order l/\q\ and acquires Coulomb divergence. In the case of the vector 
counterparts, the leading contributions to luiH and ImQ are both Coulomb divergent. 



6. NRQCD corrections 



In this section, we compute the NRQCD amplitudes of order a^. As shown in section Isl the 

andi^^) 

Cl OCl 

in the IR regions. In order to determine the short distance coefficients pj^^ and Vn^^ in 



one-loop QCD corrections "i-^^ i^ and ^-4.^^ contain Coulomb and logarithmic divergences 



(|2.7| ), which are insensitive to the long distance interactions, we remove those divergences 
based on the fact that NRQCD amplitudes must reproduce the corresponding full QCD 
amplitudes in the IR regions because NRQCD is a low energy effective field theory of 
QCD. We shall find that the divergences of iA'^J'^^^^ and i-^^^^ are identified as the one-loop 
corrections to the perturbative NRQCD matrix elements {0\On\bci)^^^ and (0|O^|5ci)^^^ in 

dU). 

Instead of following the direct NRQCD approach, we compute NRQCD quantities 
from the full QCD expressions iA'^-^^^ and iAj^^^ based on the method in |3^]. First we 
carry out the integration over the temporal component of the loop integral and, then, 
expand the integrand in powers of q/mq and k/mq, where k is the spatial component of 
the loop momentum. We regularize divergent integrals dimensionally and drop scaleless 
power-divergent integrals. The only nonvanishing divergent contributions are, then, either 
logarithmic or Coulomb divergent. As in [31|, we use a special notation for this prescription 
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for the loop integration as -J^. Once integral has been evaluated by contour integration, 
then the remaining integral is denoted as 

where indicates that the integrand of the spatial loop variable must be expanded in 
powers of q/rnq and k/ruQ and then regulated dimensionally in d — 1 spatial dimensions. 

To evaluate the vertex corrections j^^j^q^j-, and A^j^q^^j-, in NRQCD, we begin with 
the full QCD expressions in (|5.4|) by replacing the loop integrals J^. with Jfj^: 

2^ J'rV{p-b)'^^{k,P,q,mb,-mch5u{pc) 



A5,NRQCD = -^asCpf^ , (6.2a) 



Anrqcd - -'9s^^% D^D;D2 ' ^^-^^^ 



where T^{k,p,q,mb,m^ is defined in ( ^.51 ). In the NRQCD case, we omit the tensor 
reduction and directly evaluate the integrals because the tensor reduction does not 
simplify the intermediate steps of calculation considerably. Then the expressions in ( |6.2| ) 
becomes 

A^,NRQCD = -i9lCFv{Pl)\[{d - 2)5i - 2(V - < - m^S^ + 8q,S^,]j'' 

+2(7715707, - 777,7,70)^3^ + 2(2 - (i)50^7. - 89053^7, |757/(p.), (6.3a) 
AkRQCD = -ig!CFv{p-b)!^[id - 2)Si - 2(4/ -ml- ml)S2 + 8q,S',]Y 

+2{mbYlu + m,j,f)S^ + 2(2 - d)5f 7, - 8g^53^7.}^^(Pc), (6.3b) 
where the loop integrals Si, S2, Si^, and Sj^"^ are defined by 

Next we carry out the integrals in ( |6.3D by contour integration. In order to identify 
the residues of the k^ integral, we express the denominator factors of the integrands in 
( |6.4| ) in the following form: 

Do = {ky - fc2 + = (feO - |fc| + ie){k^ + |fe| - ie), (6.5a) 
Di = {k^ - Ebf -Al + ie = (A;° + Ab - Eb - ie){k^ -Ab-Eb + ie), (6.5b) 
D2 = {k° + Ecf -Al + ie = (fcO + Ac + Ec- ie){k° - A, + + ie), (6.5c) 
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where Ag for Q = c oi b are defined by 



A. 



+ {k + qf. 



(6.6) 



The resultant integrands of the {d — l)-dimensional integrals Jfj^ are then expanded in 
powers of q/raQ and k/mq for Q = b or c. The following relations are useful in these 
expansions: 



A, 



Er 



A - A 



A;2 + 2fe • q 

2 



ml 



a2 



{EQ±\k\f = ^2\k\{EQ^q-k). 



(6.7a) 

(6.7b) 
(6.7c) 



It is evident from ( |6^ ) that the factors such as 1/(Aq - Eq) and l/[Ag - (E'g ± |fc|)^] 
may give rise to IR singularities. In this step, scaleless integrals that are power divergent 
in the UV regions are neglected under dimensional regularization. The nonvanishing ele- 
mentary integrals no, rii, 77-2, and 713 that survive in this step are evaluated in appendix^. 
Nonvanishing scaleless integrals are logarithmically divergent, which are proportional to uq 
defined in ( |C.2D . The integrals ni, 712, and 77,3 defined in (|C.3| ) have scale dependencies on 
\q\. Eventually, all of the loop integrals in ( |6.4D are decomposed into linear combinations of 
these elementary integrals. The resultant values for the integrals (|6.4| ) are given in ( |C.1C| ), 
( C.16| ), ( p. 20 ), and ( p. 23 ) of appendix]^. Substituting these values to (|6.3D , we find the 
multiplicative vertex correction factors in NRQCD as 



A 



5,NRQCD 



UsCf Pi ■ Pc 
4lT 2p2 



A 



NRQCD 



+ log 



)Li( 

2{mb + rricf 



eiR euv/ d 6 



+ 6 + 



[A 



sJnrqcdI 



Pb'Pc- mbrricj 2pi ■ p, 



1 

eiR 



(6.8a) 
(6.8b) 



Substituting the multiplicative vertex correction factors in ( |6.8D and the heavy quark wave- 
function renormalization factor in (4.12) into ( [4.11 ), we obtain Gs^nrqcd and Gnrqcd as 



G 



5,NRQCD 



1 + 



OsCf pi ■ Pc 



G 



NRQCD 



2p2 

T 

= Gs^NRQCDlm, 



1 

eiR 

— +log 

eiR V 9 



+ 6 + 



4p2 
Ph-Vc 



+ 



Li(5, Cc) - 
2(r7T,{, + m-c 



Pl-Pc- mbmcj 2pi-p, 



vr 



(6.9a) 
(6.9b) 



As was expected, the logarithmic and Coulomb divergences in the IR regions of G5 and G 
in ( ^.11 ) are reproduced in Gs^nrqcd and Gnrqcd in (6^), respectively. We notice that 
unlike G5 and G, Gs^nrqcd and Gnrqcd contain logarithmic UV divergences. In a similar 
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way, the remaining NRQCD correction factors are obtained as 



Hi 



5,NRQCD 



m5 



2(mfe 



Q 



5,NRQCD 



-ffNRQCD 

Qnrqcd 



rrir 



"^{Pb ■Pc + mbrric) 



Pb-Pc- marine 
-f^5,NRQCDL^^_, 
Q5,NRQCdL,^^_, 



.10a) 



(6.10b) 

(6.10c) 
(6.10d) 



All of the quantities in ( |6.1C| ) are pure imaginary. According to (^^) and (^^), ^^5,nrqcd 
is of order \q\ and, therefore, finite and the leading contribution to Qs^nrqcd is of order 
l/|g| in the limit \q\ — )• and acquires Coulomb divergence. In the case of the vector coun- 
terparts, the leading contributions to -ffNRQCD and Qnrqcd are both Coulomb divergent, 
which show the behavior oc l/|q| in the limit |g| — ?■ 0. The expressions in ( |6.10| ) reproduce 
the IR behaviors of the full QCD counterparts in ( |5.12 ). 



7. Results for the short distance coefficients 

In this section, we list our final results for the short distance coefficients in , On , and bn 
for J = and 1 and for n = 0, 1, and 2. We have shown that the IR behaviors of G5, H^, 
Q5, G, H, and Q in ( 5.11 ) and ( 5.12| ) are exactly reproduced by the NRQCD counterparts 
Gs.NRQCD, -ff5,NRQCD, Qs.NRQCD, Gnrqcd, -f^NRQCD, and Qnrqcd in (|o|) and (|6.1C| ), 
respectively. Therefore, all of the quantities AG5, AQs, AG, AH, and AQ in ( 4.1^ ) 

are free of IR singularities. Because the imaginary parts of the QCD amplitudes are the 
same as those of the NRQCD counterparts, all of the quantities AG5, AH^, AQ^, AG, 
AH, and AQ in ( [4.18 ) are real. Except for AG5 and AG, which have logarithmic UV 
divergences originated from the NRQCD factors Gs^nrqcd and Gnrqcd, all of the other 
quantities {AH^, AQ5, AH, and AQ) are finite in both UV and IR regions. We renormalize 
the UV divergences of AG5 and AG according to the MS scheme. Our final results for 



AG 



(1^ 

5,MS' 



AH 



(1) 



AG 



(1^ 

5,MS 



and AQg 

asGp 



(1) 



are 



+ 



An 

12(6, Cc 



4 + 2 



3 + 



2 

log— I 
mf 



+ 



{nib + ruc)'^ "I 
Pb-Pc- mbrric 
Pb-P' 



Pb ■ Pc f. 



2p^ V 



2K{6,ec) 



4p2 



Li{6, Cc) - 1 



log- 



At 



mt 



mrmf. 



log- 



m 



AH, 



(1) _ asCp 



27r 
An 



+ 



rub 



Pl-Pc + mbrric 
4 



(5^Li((5, Cc 



mb + rric, m, 
+ — ^^5 — log — 



2 1 



{ruh - rric) 



mt 



b 

(7.1a) 
(7.1b) 



Pb-Pc- 

rub + mc 
2^ 



rribmc 
'^"1 ruf, 



2p^{pi ■ Pc + mbmc 



+ 



{nib + mc)5'^Li{6, Cc 



rrir 



Ap^ 



1 "T-c 
log— I 



rub + mc 
p2 



(7.1c) 
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The results for the vector part can be obtained by replacing nic in (7J) with —rric as 





= AG^^^I 


MS 


5. MS Imc^-nic 








lr?ic— >— rrtc 


AQ(i) 


= AO^^^I 



(7.2a) 
(7.2b) 
(7.2c) 



Taking mi, = rric = mg in (7.2) for the vector current, we recover the corresponding results 
of the electromagnetic current for the spin-triplet S'-wave QQ pair of the same flavor in 

Now we are ready to obtain the short distance coefficients Pn , o-n , and bn for n = 0, 
1, and 2 and for j = and 1. The order-a^ short distance coefficients can be found from 
the expansion formulas in ( [4.16| ) and ( [4.17 ) as 



Pr 



(0) 



P, 



(0) 



p 



(0) 



,{o) 



,(0) 



,(0) 



do) 



l(0) 



1, 



(mfo + me) 



-{Uml 



lOmhiric + llm^ 



1, 



(mfe 



. ^ , (llml + lOmhiric + llm 

I28mtmf. ^ ^ 



^) 



1 



2mi,mc 
4 



3(m^ + m^) 



31m^ + ISmfm^ + 31m^ 
256m^me 



(7.3a) 
(7.3b) 

(7.3c) 

(7.3d) 
(7.3e) 

(7.3f) 
(7.3g) 
(7.3h) 

(7.3i) 



It is straightforward to obtain and bn^ for n > 3 in the same way. The S- and D-wave 
components of the short distance coefficient for the vector current can be obtained from 



O) as 5, 



(i) 



aii'^ + ifti^'^ and dI^^ 



ttn' ■ The short distance coefficients of order are 



free of scale dependence. The values for Pq'^ and agree with the previous results in 



1 24, 25, ^Q*^^ and S\^' agree with the previous results in [25] 



,(0) 



The order-a^ short distance coefficients Pn^^ 
by substituting [AG^^^J^jg, AH^^\ 



and AQi^^ 



for n = 0, 1, and 2 can be obtained 



.(1) 

'5 ' 



r(l) 



and AQ^^^ in (|4lg| ), 
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respectively. The results are 



[^1 



(1)1 



IMS 



1^2 J MS 



3(m6 — rric 
2(mb + m-c 



21 



(7.4a) 



-1 



47rm^m2 j 144 



4(m^ + 98mfcmc + m^) + 3(7m5 + 46m;,mc + 7m; 



X log — £ 



mi, + rric 



2/ x2, 

+ o (™b + "^cj log 

3 mhm^ 



(7.4b) 



4(4089m^ - 10364mbmc - 36586mf - 10364mbm^ 



QsCf [ -1 
47rm^m4 \ 57600 

+4089m^) - 15(177m^ + 2212m||mc + i582mlml + 2212mbm^ + 177m^) 



rub — rric , 

X log — £ 

m 



mb + rric 



b-i 



{nib + mc 
60 



-{21m1 + 2mbmc + 21m^) log ■ 



.(7.4c) 



Except for Pq^^ , P,!^^ for n > 1 are dependent on the NRQCD factorization scale fi that has 
been introduced in the process of renormalization. Our result for Pq^^ agrees with those in 



13, H]. The results for [P| 



(1)1 



iMS 



and [-fi^'^lMs new. 



The short distance coefficients and bn^ for n = 0, 1, and 2 of the vector current 
at order a] are obtained by substituting [AG(^)]j;jg and Aij(i) into AG(i) and AHW 

m 



(4.20), respectively. The results for ai^^ are 



(1) _ asCp 

OIsCf 



n Jms 



■*2 Jms 



3(mfc — rric 
2{nib + rric 



log' 



rrr" 



rrrf 



1 



47rm^m2 [ 144 



32(m^ — mbrric + m^) — 3(7m,^ + lOmfcrn-c + 7ml 



mb-mc 
mb + rric 



21 



log -4 



rrrf 



Airmlmi \ 19200 



(7.5a) 



2 1 

+ -(rrrb + mc)^log^— L (7.5b) 
o mbmc J 

16(547m^ + 128rrrf r^c + 122mlml + USmbm^ + 547rrr^) 

2 



-15(59m^ + IMmlmc + IMmlml + 164rrrbm^ + 59m^)^^^^— ^ log 



mb + m-c 



1 

■ — {mb + mcf{7ml - Qmbm^ + 7rrr^)log 

20 mbmc 



(7.5c) 
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And the short distance coefficients bn^ are 



y"i J MS 



[4 ^Ims 



UsCf / -1 
A-Kmlml \ 288 



2+ ."^^""^Mog^ 



4(mb + TTic 



8(19m2 - lOm^m, + IQm^) - 21 (mf + 4m6m, 



ruh — rric , 
X log — 

ml 



-(mfe + mc) log 



mhUic 



47rm5m5 1 115200 



1 



(7.6a) 



(7.6b) 



8(5647m^ - 972mgmc + 922m^m^ - 972mbml + 5647m^) 



2n 



-15(527m^ + 1912m^mc + 2082m^m2 + 1912mfcm^ + 527^^)^^^^-^^ log ^ 



mi, + nic 



mt 



120 mfjiTic 



(7.6c) 



As in the case of and ft^'^^ , it is straightforward to obtain and bn^ for n > 3. Except 
for Cq^^ and 6q^\ ai^^ and fti^'* for n> 1 are dependent on the NRQCD factorization scale 



On^^ + ^b^Q ^ agrees with that in 



The results for [a-[ 



(1)1 



IMS' 



jj,. Our result for — u,q 

[^1^]mS^ [4^^] MS' [4^^] MS 

Let us make a rough estimate of the effect of relativistic corrections at order . It is 
convenient to define the ratio of the NRQCD matrix element of relative order-g^" to the 
leading-order matrix element: 



{0\On\Bc 

{0\Oo\B, 



{0\0\o\B*c 



(7.7) 



where On and are defined in (|4.1|), and we have used the property that {q )b* is 
independent of i. The ratios {q^^)Bc,B* are normalized to be consistent with that, (q^")^, 
for the QQ quarkonium H considered in [39, 40, 31]. In p^, it was shown that the MS 



value for {q'^^)H satisfies a generalized Gremm-Kapustin relation |41|: 



lie)] 



MS' 



(7.8) 



We assume that this relation is still valid in the case of the be meson. Unfortunately, 
unlike the case of the 5-wave bound states of the QQ pairs, the leptonic decay rates for 
Be and B* have not been measured accurately so that one cannot determine {q^'^)Bc,B* 
with empirical data. Instead, by taking upper and lower bounds of the ratios {q^^)B^^B* 
as {q^)j/^ = 0.441 GeV^ |o| and {q^)r(is) = -0.193 GeV^ ||, respectively, we make a 
rough estimate of the sums of products of 5-wave short distance coefficients and operator 
matrix elements. Taking the central values for the variables mf, = 4.6 GeV, rric = 1.5 GeV, 
// = (mb + m,)/2 = 3.05 GeV, and {q^)B.,Bt = {{q^)j/^ + (g')T(i5))/2 = 0.124 GeV^, we 
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find that 
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n=0 
1 

E 

n=0 
2 

E 

n=0 

oo 

E 

n=0 


E 

n=0 
1 

E 

n=0 
2 

E 

11=0 
oo 

E 

n=0 



rp(i)" 

^ n 


MS ^ 


rp{i)" 


MS ^ 


n 


MS ^ 


^ n 


MS ^ 




[ 

MS 




[ 

MS 




[ 

MS 


[9(1)" 


[ 

MS 



(9^)bc]Sis 



/Be 



MS 



asCp 
47r 



X 4.292, 



X 4.293 



-0.185 
-0.330' 



X 4.294+0-^52 



0.306' 



MS 



X 4 294+°-^52 



OisCp 



X 6.209, 



X D.±DO_Q 289, 



X 6.170;[J;2i3, 



IMS 



V fi 170+0.213 

4^ XD.im_Q 252, 



(7.9a) 
(7.9b) 
(7.9c) 
(7.9d) 
(7.9e) 
(7.9f) 
(7.9g) 
(7.9h) 



where the uncertainties are estimated by varying the values rric < < irih and —0.193 GeV^ 
— (q^)bc ^ 0.441 GeV-^. The estimates in ( |7.9D show that the series expansions for the 
relativistic corrections to the axial vector and vector currents at order as converge rapidly in 
spite of the large uncertainties in the ratio {0^)3^- addition, the order-one contributions 
of the form (mfe — mc)/{mi, + mc) log(m^/m^) do not deteriorate the convergence of the 
short distance coefficients for the axial vector and vector currents. In |26], the authors 
obtained the short distance coefficient of order-a^w" as —24.4 x [as(m;,)/7r]^ with the 
input parameters m^y = /i = 4.8 GeV and ttt-c = 1.65 GeV. If we use these values for m^, 
TTic, 11 and set {q^)Bc = 0.124 GeV^, then the coefficients ( [7^ ) vary by about < 3%. As an 
example, ( 7.9d| ) becomes 



E [pi" 

n=0 



MS 



MS 



4tt 



X 4.399. 



(7.10) 



8. Summary 

We have computed the short distance coefficients for the NRQCD factorization formu- 
las of the meson-to-vacuum matrix elements for the axial vector and vector parts of the 
charged weak current in the 5"- wave spin-singlet and -triplet be mesons, respectively. The 
computation was carried out at order and a] including relativistic corrections of all of 
the be NRQCD operators that contain any number of ordinary derivatives without gauge 
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fields. In the Coulomb gauge, gauge field contributions first appear from relative order 
Vq. The numerical value for the correction of order agVQ is tiny (<0.1% of the leading 
order contribution). We have reproduced all available short distance coefficients of order 
as or less and our results of order al with relativistic corrections are new. By taking the 
limit nih = rric = rnq, we have reproduced the results for the order-a^ corrections to the 
quarkonium electromagnetic current for the spin-triplet S'-wave QQ pair with the same fla- 
vor. Although we have listed explicitly the first few short distance coefficients of order-a^ 
for the relativistic corrections, it is straightforward to obtain the terms of higher orders 
in g^". The results reveal that, in spite of the large uncertainties in the ratios {q^)Bc,Bi-, 
the relativistic corrections to the leptonic decay rate of the Bc{Bc*) meson at order a] or 
less converge rapidly, which shows a strong contrast to the uncomfortably large corrections 
of order 0^191*^. The short distance coefficients listed in this paper can be employed to 
compute the resummation of relativistic corrections to the phenomenological measurables 
that involve Be and B* production and decay. 
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A. Tensor-integral reduction 

In this appendix, we describe the tensor-integral reduction that we use to simplify ( ^.4| ) to 
obtain (|5.6|) . 

Tensor integrals of rank 1 and 2 that depend on p or on both p and q can be expressed 
in terms of scalar integrals as follows: 

k'^fik^p) = ^ [ p-kfik,p), (A.l) 
k P Jk 



k^'k''f{k,p) = g^^ / di{k,p)f{k,p) +p^'p- / d2{k,p)f{k,p), (A.2) 

Jk Jk 

k''f{k,p,q) =p^ [ d3{k,p,q)fik,p,q)+qf' [ d^{k,p,q)f{k,p,q), (A.3) 
Jk Jk 

k^^k'^f{Kp,q)=g^^ [ d^[k,p,q)f{k,p,q)+p^p^ [ deik,p,q)f{k,p,q) 
Jk Jk 

+q''q'' / d7{k,p,q)f{k,p,q) 
Jk 

+ {p^q-+p-q^^) [ dsik,p,q)fik,p,q), (A.4) 
Jk 
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where k is the loop momentum, the symbol is defined in ( |5.2[) , and / is an arbitrary 
scalar function of the argument four- vectors. The functions dj's are defined by 



di{k,p) 
d2{k,p) 
d^{k,p,q) 



db{k,p,q) 
d(>{k,p,q) 



1 



d-1 
1 



{k-pf 



k ■ qp ■ q — k ■ p 



(p ■ qY — p'^q^ 
d4{k,p,q) = d-i{k,q,p), 
1 



d-2 



^2 _^ /(g • fc) ^ + -kf -2p- qp-kq - k 



{p ■ qY — p^q^ 



(A.5) 

(A.6) 

(A.7) 
(A.8) 
(A.9) 



1 



[(A; • pf{d - 1) - p^k^]q^ + {d-2){p- qf{k ■ qf 



{d-2) [{p-qY-p\^f 
+ [k\p ■ qf - 2{d - l){k ■ p){k ■ q){p ■ q) + p\k ■ qf]q^^, (A.IO) 
d7ik,p,q) = d6{k,q,p), (A. 11) 

ds{k,p,q) = , , \, [q^{k ■ pf - k'^ip ■ qf + p\k ■ qf + p'^q^k'^] p ■ q 



{d-2)[{p.qf-p^q^f 
2p^q^ k ■ pk ■ q + d\p ■ qk ■ p — p^ k ■ q\ \p ■ qk ■ q — q^ k ■ p\ |. (A. 12) 



B. Scalar integrals for the vertex corrections 



In this appendix, we list the definitions and the values for the scalar integrals Jj that appear 



in the vertex corrections in (5.6). The scalar integrals Jj's are defined by 



J,; 



(B.l) 



where k is the loop momentum, the symbol and the denominator factors DiS are defined 
in (|5.2| ), and the numerators A'j's of the integrand are defined by 



N5 



1, 

2k -q, 
1 

d-2 
2k -p, 



^2 ^ • kf + q'^jp ■ kf - 2{p ■ q){p ■ k){q ■ k) 
[p ■ qf — p^q^ 



(B.2) 
(B.3) 
(B.4) 

(B.5) 

(B.6) 
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A^6 



+ 



nih + rric q 
d — 2 (p ■ — p^q^ 

rUc - m-b 



2 , {d-l)\p-'{q-k)-{p-q){p-k)f 



{p ■ qY — P^q^ 



rric + mh 



-k^{p-q) + {p-k){q-k) + 



x\p^{q-k)-{p-q){p-k)] 



{d-l)[q\p-k)-{p-q){q-k)] 



{p ■ ~ p'^Q^ 



Ny = 
+ 
+ 



nib + w-c p'^q^ 



k\p ■ q) + (p ■ k)iq ■ k) 



d — 2 (p ■ qY — p'^q^ . 
{d - l)[q^{p ■k)-{p- q){q ■ k)]\p'^{q ■ k) - {p ■ q){p ■ k)] 



rric - mb 
rric + irib 



k'q' - {q ■ ky + 



{p ■ ~ p^Q^ 

{d-l)[q^{p-k)-{p.q){q-k)f 



(p ■ qY — p^"^ 



(B.7) 



(B.8) 



where d = 4 — 2e is the number of space-time dimensions. The external momenta p and q 
are defined by p = \ij)c a-nd q = \{j>c — Pi)-, where p^ and Pc are the momenta for the 
h and c, respectively, which are on their mass shells: p| = and = fr?c- 

The values for the scalar integrals Jj's are 



Ji 



(47r)=^ 



euv 



+ 2 + 2mS - 25^Li(5, Cc) + ec log 



+(l-ec)log 



47r/i^e 
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(47r)2 8p2 V eiR 2 



y - -Li((5,ec) 
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+ 2if(,5,ec)-y 
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+ n-^2('5, ec)l0g 



2 I ' 



J4 = 



= 



(47r) 



(47r)^ 



2z7r(5- 2(5^Li((5, ec) + 



+ Ji 



8p2 



2(47r) 



log 



2 ' 



i (rrife + mc)q^ 



(47r)2 pi-pc- mbiric 



S^MS, ec) -inS+\( - <^ ] log ^ 



4 V mb + mc 



4p2 



(B.9) 



(B.IO) 
(B.ll) 
(B.12) 
(B.13) 
(B.14) 
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(mb + TTlc 



4(47r)2 pi' Pc — rrihrn. 
+ 



ml-ml + 2{pi ■ pc - mbfric) 



rub - rric 



nib + rric 



ml — ml / nib — nT-c "f^l ~ "^c 



mb + mc 



4p2 



Pb- Pc- rubmc 



mi — m. 



2' 



ITT' 



4p2 



xlog— I }, 
ml I 



(B.15) 



where the variables 5 and Cc are defined in ( |3.4D and the functions Li((5, Cc), L2{5,ec), and 
K{5,ec) are defined in (|5.9D. 



C. Integrals for NRQCD vertex corrections 

In this appendix, we hst elementary loop integrals that are useful in computing the NRQCD 
corrections considered in section ^. We also list the values for the integrals defined in (^^) . 
We follow the strategy of evaluating the integrals given in |31]. 

C.l Elementary scalar integrals 

In dimensional regularization, scaleless power-divergent integrals vanish: 

1 







for n 7^ 3. The only nonvanishing scaleless integral is 



no 



1 



1^1 



1 

4^ 



1 



euv eiR 



which diverges logarithmically. Except for the integral (| 
depending on \q\, which are 

1 i 

ni = 



1^ k"^ + 2k ■ q — ie An 
1 



(C.l) 
(C.2) 

nonvanishing integrals are 
(C.3a) 



n2 



fc k^{k^ + 2k-q-ie) 



1 



ns = 



167r|g| VeiR 

fc2 



h log 



(C.3b) 



(C.3c) 



1 1^ k"^ + 2k ■ q — ie 2tt 
In projecting out the S'-wave contribution from a scalar integral that depends on q, 
we have to take the average over the angle of q. The following formulas are useful in that 
step: 



f{k' 



1 



kEQ±q-k 2|q| 



log 



1 



{Eb-q-k){E^ + q-k) Eb + E^Jk 

1 



2{Eb + Ec)\q 



Eq + \q\ 
Eq - \q\ 

log 



f{k'), 
1 



(C.4) 



1 



+ ■ 

Eb + q ■ k Ec — q ■ k 

iEb + \q\)iE, + \q\) 
(Eb - \q\)iE, - \q\) 



[ f{k% (C.5) 

Jk 
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where Q = b ov c, k = k/\k\, and f{k'^) is any function of fc^. 

In the fohowing sections, we express the integrals Si, S2, S!t, and S^'^ defined in 



in hnear combinations of no, ni, 712, and 113 in ( C.2 ) and ( p.3D . We also find the covariant 
forms of the integrals Si, S2, Si^, and S^"^ . 

C.2 Si 

The integral of Si defined in (^^) is the sum of two contributions: Si = Sic + Sii, 
where Sic and 6*^5 are the contributions from the poles of the charm quark and antibottom 
quark, respectively. 

The contribution from the c pole is 



'Ic 



(C.6) 



We find that all of the factors in the denominator of the integrands are of order mg as 
q — )• and — )• 0. Therefore, the expansion of Aq in powers of {k + q)^ /rnq gives only 
scaleless, power-divergent integrals, which vanish, so that 

Sic = 0. (C.7) 

The contribution from the h pole is 

m + Ec)%k^ + 2k.q-ie- ^^'^^ 
Expanding Aq in powers of {k + q)'^ /tjiq, we obtain 

By adding the two contributions (iC.Tl) and (|C.9| ), we obtain 

i 



where 5 is defined in 



C.3 S2 

The integral of S2 is the sum of three contributions: S2 = S2g + S2c + S2i, where S2g, S2C, 
and S2I are the contributions from the poles of the gluon, charm quark, and antibottom 
quark, respectively. 

The gluon pole contribution is 



S2g 



8Jk \k\^{Eb-q-k){Ec + q-k) 

{Eb+\q\){Ec+\q\ 



16{Eb + Ec)\q\ 



Uq log 



{Eb-\q\){Ec-\q\ 



(C.ll) 
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which is proportional to the scaleless logarithmicahy divergent integral uq. 
The contribution from the c pole is 



S' 



2c 



1 



4iEh + E,)Jk A,(Ae + EcW - (Ac + E,)^ - ie] ' 



(C.12) 



We find that the factors Ac and Ac + Ec are of order niQ as g — )• and fc — )• 0. Therefore, 
the expansion of the factors 1/Ac and l/(Ac + Ec) are trivial and the expansion gives only 
scaleless power-divergent integrals. The expansion of the last factor — {Ac + Ec)"^ — ie] 
can be done in powers of fc^/(Ac + -Ec)^ and then the factor l/(Ac + Ec)"^^ is expanded in 
powers of {k + qf' /mq. We find that all of the contributions are scaleless, power-divergent 
integrals so that 

S2c = 0. (C.13) 



The contribution from the h is 



S. 



2b 



1 + {Eb/Ab 



A{Eb + Ec) h [k^ - ( A6 - Eb)^ - le] [k^ + 2k ■ q - ie] ' 



(C.14) 



The expansion of the integrand for is similar to that used in deriving (67) of pi| . 
Following that method, we find that 



S. 



2b 



-71-2 



1 ^ vr^^e 

+ log T\ h ivr 



2{Ei, + Ec) 2,2TT{Eb + Ec)\q\ \eiK q 

The sum of the three contributions in ( p.ll| ), ( p.l3| ), and ( |C.15 ) is 



(47r)2 8p2 



euv eiR 



Li{6, Cc 



vr it: 

T ^ T 



1 , vr/i^e-'^E 
— + log 
eiR q ^ 



(C.15) 



(C.16) 



where 5 and Cc are defined in ( ^^ ) and the function Li(5, ec) is defined in (^). 
C.4 

The integral is the sum of three contributions: = Si^^ + S^^ + S^^, where S^^, 



and 5*2^ are the contributions from the poles of the gluon, charm quark, and antibottom 
quark, respectively. 

Following the same way that has been used to evaluate the integrals of 5*1 and ^2, 
we carry out the integrals for S^^, 5*3^, and by contour integration. The sum of the 
three contributions is 



cO 
-^3 



8/fc \k\^{Ei, - q ■ k){Ec + q-k) 4(^6 + Ec)Jk 
1 



Ac[A:2 - (Ac + Ec 



Afc[fc2 _ (Afc - Eh)^^ 



ie\ 



le] 

(C.17) 



SI 



k^ 



8/fc \k\3{Eb-q-k){Ec + q-k) ^E^ + Ec) Jk 



k' 



Ac(Ac + Ec) 



[1 + {Eb/Ab)]k' 



fc2 _ (Ac + Sc)2 - ie 



[fe2 - (Afc - Efc)2 - ie] [fc2 + 2fc • q - ie] 



(C.18) 
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^3 and SI include three terms, which correspond to ^g^, S'g^, and Si^-^, respectively. 

The integrands of the gluon pole contributions 5*3^ and S^g, which are the first terms 
in 5*3 and S^, respectively, have factors 1/{Eq ± q ■ k) that expand in powers of q ■ k/Eq 
producing scaleless factors and the factor does not generate logarithmic divergence. 

Therefore, Si^g = 0. The second terms of 5*3 and S"! are the charm quark pole contributions 
^3^ and S\^. We can follow the same procedure that was employed in expanding the 
integrand for 820 in ( |C.12 ) to find that S'^^ = 0. The last terms in 5*3 and S\ are the 
antibottom quark pole contributions and 5*^, whose structure is similar to that of S2I 
in ( C.14 ). We find that S^j^ contains only scaleless, power-divergent integrals, which vanish, 
and the only nonvanishing contribution is 5*^: 



^36 



■ni 



4{Eb + E,)q^'^' W7r{Eb + E^ 
The Lorentz covariant expression for Si^ is, then, obtained as 

p-q 



ITT 



(47r)2 4(5p2 



-p^ + 



(C.19) 



(C.20) 



C.5 ST 



The integral S'^^ is the sum of three contributions: Sf^ = S^g + S^^ + S^^, where S^^ , S^^ , 



46 



■'Ag 



and S^^ are the contributions from the poles of the gluon, charm quark, and antibottom 
quark, respectively. 

After evaluating the integral by contour integration, we find that 



si' 
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k \k\{Eb-q-k){Er^ + q-k) 
1 - (Eb/Ab) 
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A{Eh + Ec)Jk 



1 + (EJA,) 
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8Jk \k\^{Eb-q-k){Ec + q-k) 
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A{Eb + E,)fk 



fc2 _ (Ae + Ecf 



l£ 



(C.21a) 



{Ak-E,f 

k^y 



ie\ 



k \kmEf,-q-k){E, + q-k) 
[1 + {E,/A,)]k'k^ 



+ 



A{Eb + E,)Jk 



A,[k^-iAc + E,)^-ie] 

(C.21b) 

k'k^/[A,{A, + E,)] 



[fc2 _ + E,)^ 



l£\ 



[fc2 - (Ab - EbY - ieW + 2k q 



(C.21c) 



Like S'g, the three terms in each of S^^, S^^, and S^^ correspond to 5*4^, S^^ , and S'^^ , 
respectively. Following the same procedure to expand the integrands of S'g , we find that 



2{Eh + E,){d - 2) 



32Tr{Eb + Ec 



5'^ + 



q'q^ 



ni 



4g2 



rii 



{d-l)n3 
4q2 



(C.22) 
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The Lorentz covariant expression for S^'^ is now obtained as 



5r 



1 



452^4 



+ 



p2 452^2 4j2p- 



(47r)2 2 
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